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Abstract 

Many observables of interest in lattice QCD are extracted from correlation functions involving 
the vector current. If Wilson fermions are used, it is therefore of practical importance that, besides 
the action, the current be 0(a) improved in order to remove the leading discretization errors 
from the observables. Here we introduce and apply a new method to determine the improvement 
coefficient for the two most widely used discretizations of the current. 
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I. INTRODUCTION 


Lattice QCD is a powerful tool to calculate the predictions of Quantum Chromodynamics 
in its non-perturbative regime. While the quantum held theory is regularized by discretizing 
it on a lattice, ultimately the quantities of interest - for instance, ratios of hadron masses - 
must be determined in the limit where the cutoff is removed. For numerical purposes, it is 
computationally advantageous to accelerate the approach to the continuum by removing the 
leading cutoff effects. In particular, Symanzik’s continuum effective theory BE] can be used 
to remove the 0(a)-cutoff effects which appear generically when using the Wilson fermion 
action in lattice QCD simulations [3]. To eliminate 0(a)-cutoff effects in the hadronic 
spectrum it suffices to improve the action by introducing the dimension-hve Sheikholeslami- 
Wohlert term |lj with a non-perturbatively determined coefficient Csw [S]. However, the 
addition of higher-dimensional counterterms to local operators is also necessary for the 
improvement of their matrix elements, along with the determination of the corresponding 
improvement coefficients. 


In the following, we focus on the vector current, which requires a single 0(a)-improvement 
term. Estimates of the relative contribution of the improvement term evaluated with the 
perturbative improvement coefficient may suggest that the effect of the improvement in cor¬ 
relation functions would be small for the local vector current [6]. However, an improvement 
condition based on chiral Ward identities previously used to determine the improvement 
coefficient Cy non-perturbatively in pure gauge theory with |6] and without Schrodinger 
functional boundary conditions indicated signihcant deviations from the tree-level result. 


In this work, we describe a simple prescription for the non-perturbative determination of 
the improvement coefficients, Cy and Cy, for the local and conserved (point-split) isovector 
vector currents, dehned below. In the following section we report large differences between 
the lowest-order perturbative estimates and our non-perturbative evaluation of the improve¬ 


ment coefficients with W = 2 Wilson clover fermions. In section |IV| we demonstrate the 
effects of the improvement on the scaling of an observable in the continuum limit. 
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II. THEORY BACKGROUND AND A NEW IMPROVEMENT CONDITION 

We use the 0(a)-improved Wilson fermion action with the non-perturbatively tuned value 
of Csw [B]- The two discretizations of the continuum vector current that we employ are 


^ {^{x + a/i)(l + li,)Ul{xf-^i){x) 

-^(a;)(l - l^,)U^{x)^^p{x)^ . (2) 

The renormalized improved current for i = /, c is dehned by [5] 

(V'r);(i) = Zlj(l+ tVom,)(V'i);(i), (3) 

('''i);(i^) = (r)i.(x)+acVS,U„(i), (4) 

where the lattice discretization of d^T^y{x) will be discussed later and in the 


notation of [5], while ^^ = 1. The on-shell improvement of the vector current is required 
in many lattice studies: hadronic form factors, the hadronic contributions to {g — 2)^ and 
thermal correlation functions related to the dilepton production rate, to mention a few. 

While the local vector current requires improvement only at one-loop order [H] the con¬ 
served vector current requires improvement at tree level also in the massless limit. 


0.01225(1) X Cf X 4 + 0(4), 

(5) 

+ 0(4), 

(6) 


where Cp = — 1)/(2V) is the quadratic Casimir in the fundamental representation for 

gauge group SU(V). 

A. Improvement condition for Cy^ 

The proposed improvement condition is based on two discretizations of the vector current 
dehned in eq. Q (/) and eq. (|^ (c). The main observable we consider is the vector current 
correlator 

{Gi)U^o) ^ {Gv)%{xo) 

+ Z^Z^ j d3x(a4(U);(xo,aj)a,T,,(0) 

+ ac\,dpT,p{xo,x){VyMy (7) 
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FIG. 1. Unimproved correlators for ensemble A5 (left) and N6 (right). 


with 


(G„)« (^o) = 44/d“^((4:,4,ii)(r)j(o)). (8) 

where we have indicated its dependence on the discretization of the current at finite lattice 
spacing, ij = ll,d,cc. We will use periodic boundary conditions in space and thermal 
boundary conditions in time, though our method is more generally applicable, for instance 
to open boundary conditions in time. In ref. [10], two discretizations of the unimproved 
vector current correlator demonstrated signihcant differences in the region of f ^ O.Sfm at 
intermediate lattice spacings corresponding to bare lattice couplings /3 = 5.3 with N{ = 2 
Wilson clover fermions. This suggests that the three independent discretizations of the 
temporal vector current correlator could be used to formulate an improvement condition. 

In figure]^ we illustrate the discrepancies between the three discretizations for two lattice 
spacings corresponding to bare lattice couplings of /3 = 5.2 and /9 = 5.5 in the left and right 
panels respectively. The details of the ensembles and the number of measurements are 
listed in table The non-perturbative renormalization constant, Zy, is taken from ref. HU. 
Note that a more precise non-perturbative result for has been reported in ref. |T2]. As 
expected, the differences between the two discretizations are reduced as the lattice spacing 
decreases. 

Demanding that three discretizations of the temporal vector current correlator agree at 
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a particular Euclidean time xq, 


(Gi)'U^o) = (GjrJ.M = (Gj)ZM, (9) 

allows one to solve the following 2x2 linear system for the improvement coefficients Cy‘^(xo)'- 

2ZyGgixo) - Gf^ixo) -Gilixo) 

-ZyGf^ixo) 2G^^(xo) - ZyGZixo) 

{GvTJ^ixo) - (Gu)L(a:o) 
iG\jYJ^{xo) - (Gu)xx(^o) 


X 


I 4 ' 

\=-i 


I a ^ 


( 10 ) 


where 


<^^( 3 ^ 0 )= / (Yx{{Vy {xo,x)dpT^p{^)). 


( 11 ) 


Translation invariance and time-reversal antisymmetry of the vector-tensor correlation func¬ 


tion was used to simplify the system (10). The correlators appearing in (10) may be averaged 


over the spatial components to improve the signal. We define the improvement coefficient 
to be ciy = cY{xq) for some choice of xq. The method is viable in practice provided that a 
signal exists both for the r.h.s. and the determinant of the linear operator on the l.h.s. of 
eq. ([^. The results for Cy obtained using different legitimate prescriptions will in general 
differ by 0(a) corrections. Ideally, one would choose Xq in a region where there is both a 
signal and higher-order lattice artifacts are highly suppressed. This improvement condition 
can be implemented directly in a hnite volume and is straightforward to compute, not re¬ 
quiring the three-point functions of ref. [6]. Although the quark mass-dependence of the 
renormalization factor is neglected in this improvement condition, namely the by term in 
eq. (|^, due to the smallness of the quark mass these effects ought to be small. Furthermore, 
in the following section numerical evidence demonstrates that no dependence on the quark 
mass is likely to be observed in these improvement coefficients. 


1. Discretization of dp and Tp^ 

In the improvement term, we use the local discretization of the tensor current, 

Tpu{x) = -^Y{x)[i,„'iY^Y{x), 


( 12 ) 
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/? 

a(fm) 

m7r(MeV) 

T X 

# cfgs 

sources 

name 

5.2 

0.079 

312 

64 X 32^ 

217 

8 

A5a 





404 

4 

A5d 

5.3 

0.063 

451 

64 X 32^ 

421 

4 

E5g 



324 

96 X 48^ 

294 

4 

F6 

5.5 

0.050 

340 

96 X 48^ 

568 

1 

N6 


TABLE I. Details of CLS Af = 2 ensembles and number of measurements used in this work. 

with the same spacetime argument as the vector current. The choice of the discretization of 

affects only higher-order lattice artifacts, which nevertheless can be large. In ref. [13], the 
improvement of the conserved current ([^ was considered. The effect of using the symmetric 
derivative du and the tensor current averaged over sites x and (x -b afi) was examined at one- 
loop level in lattice perturbation theory. While with this choice the identity = 0 still 

holds in on-shell correlation functions, it was noted to introduce large higher-order lattice 
artifacts to the connected part of the hadronic vacuum polarization tensor. Therefore, 
for the tensor current with time argument Xq (assumed positive) in correlation function 
([^, we choose the forward hnite-difference derivative, while for the tensor current at the 
origin we used the backward derivative. We remark that it is admissible to use different 
discretizations of in the determination of the improvement coefficient and subsequently 
in matrix elements of the improved operator [T3]. 

III. EVALUATION OF 

Using the improvement condition eq. ^ we determined Cy{xo) at three lattice spacings 
on ensembles with Af = 2 Wilson clover fermions with non-perturbatively tuned value of 
Csw [H] and the plaquette action. The most relevant parameters are given in table for more 
details on the ensembles, see [T3|, where the action is also given explicitly. In figure]^ we 
show the dependence of on the choice of xq for ensembles A5, F6 and N6 with bare lattice 
couplings [3 = 5.2,5.3 and 5.5, respectively, and almost identical pion masses. Additionally, 
we compare with another ensemble, E5, with (3 = 5.3 and a larger quark mass. In figure 
(left), some evidence can be seen for a plateau in the value of Cy at the smallest lattice 
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FIG. 2. Non-perturbatively determined improvement coefficient Cy (left) and Cy (right). The blue 
bursts have been displaced horizontally for clarity. 


Ensemble 

W 

AC 

Cy 

A5 

-0.434(5) 

0.203(7) 

E5 

-0.400(6) 

0.229(8) 

F6 

-0.401(6) 

0.232(9) 

N6 

-0.377(7) 

0.243(10) 


TABLE II. Results for improvement coefficients for the local and conserved vector currents for 
the ensembles listed in table The statistical error is estimated by single-elimination jackknife 
resampling. 

spacing corresponding to /3 = 5.5 (black squares). For Cy, no significant changes can be 
observed as a function of the lattice spacing. The right-hand panel of hgure shows there 
is a greater dependence on the cutoff for Cy. We are unable to distinguish any dependence 
on the quark mass. 

Our choice for the improvement coefficients is = Cy{xo/a = 3) which is used in the 
rest of this work. Our results are given in table |TT| 

This choice leads to values of Cy which deviate signihcantly from the perturbative estimate 
of eq. ([^. We can make a rough comparison with the value of Cy determined in the quenched 
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theory [6] based on the observation that, at hxed bare coupling, the improvement coefficients 
are independent of iVf to one-loop in perturbation theory. We note that the value we obtain 
at = 5.3 is quite similar to the value obtained at /? = 6.0 in the quenched theory. 

In £gure|^we show the effect of the improvement on the vector correlators with the given 
improvement condition for the A5 (left) and N6 (right) ensembles. By dehnition, the central 
values now coincide at xo/a = 3. The effect of the improvement appears to be smallest for 
the local-conserved vector current correlator, due to the contributions of the improvement 
of each current entering with opposite signs. 


2. Interpolation in Qq 

The following polynomial interpolation formulas in = 6//3 can be used to determine 
the improvement coefficients for N{ = 2 flavours of non-perturbatively improved Wilson 
fermions 

cWo) = -0.01225(1 + 7.1992 + lO.lSjo*) , (13) 

4(9o) = ^ (1 + 0 - 0.728g2) . (14) 

Conservatively, these parametrizations should be used in the interval 5.2 < (3 < 5.5, even 
though the known behavior at small is built in. It would be interesting to extend the 
present calculations to smaller values of gQ to make explicit contact with the one-loop result. 


IV. CONTINUUM LIMIT OF IMPROVED OBSERVABLES 

In order to quantify the effect of the improvement we examine the scaling of the observable 

Ii = [ dxoX^(Gi)^4(a:o) (15) 

Jtt 

and its unimproved counterpart J^, which is dehned analogously with the unimproved two- 
point function (Gu)^4(xo) (eq. ([^), toward the continuum limit. The limits ta/a = 8 
and tb/a = 26 are hxed at the smallest lattice spacing, corresponding to ta ~ 0.4fm and 
tb ~ 1.3fm. Although the contact term does not contribute to such an observable when 
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FIG. 3. Improved correlators A5 (left) and N6 (right). 


ta = 0, the lower limit explicitly removes very short-distance contributions. A lattice esti¬ 
mate for this observable is obtained from the discretized correlation function by quadrature 
with an improved integration scheme. An interpolation is needed at the limits for the coarser 
two lattice spacings. This observable is related to the slope of the Adler function through 
the mixed representation of the hadronic vacuum polarization function [16]. Therefore, it 
may serve as a useful proxy to quantify the effect of the improvement on phenomenologi¬ 
cally relevant observables which are dominated by the long-distance physics of the vector 
correlation function. 

The values of the Sommer scale, tq/q, used to set the relative scale and perform the 
continuum limit were taken from ref. na. 

In figure]^ (left) the three discretizations /((, are shown in red, green and blue, 

respectively. The scaling of the unimproved observables is modelled linearly in a to obtain 
the continuum limit. While the continuum limits of the different discretizations agree within 
the statistical precision, and could be constrained to agree in a simultaneous £t, the use of a 
single discretization demonstrates a signihcant fraction of the uncertainty in the continuum 
result is due to the long extrapolation in a. 

In the right-hand panel of hgurej^the analogous plot for the improved observable is shown. 
A quadratic model in a describes the continuum scaling well, and the error in continuum 
limit is correspondingly reduced. Furthermore, the residual scaling violations appear to be 
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FIG. 4. Linear scaling of observable (15) for unimproved currents (left) and quadratic scaling of 
improved currents (right). The points in the right-hand panel have been displaced horizontally for 
clarity. 



FIG. 5. Scaling of the ratio of local-local and conserved-conserved observable (15). 


small. Note that the observable dehned from the local-conserved vector correlation fnnction 
appears to have the mildest scaling of all three nnimproved discretizations. 

Another illnstration of the improved scaling behavionr of the 0(a)-improved observable is 
shown in £gnre[^ The ratio of the observable dehned using the local-local vector correlation 
function and the conserved-conserved vector correlation function is shown for the unimproved 
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(red) and improved currents (green). This discrepancy should vanish in the continuum 
limit and the hts (solid lines) are constrained to vanish there likewise. Here, owing to the 
additional constraint imposed in the continuum limit we can use a quadratic model for both 
unimproved and improved discretizations. 


V. CONCLUSIONS 

The basic idea of the improved strategy used here is that if discretizations of a current 
are considered, there are ndin^ + l)/2 lattice versions of its two-point function. Requiring 
the equality of these two-point functions thus allows one to determine the improvement 
coefficients of each discretization, for Ud sufficiently large. The method however requires 
that the relative normalization of the operators be known at the outset. For the program 
to go through consistently one thus needs to determine the renormalization factors in a 
situation where the improvement terms do not contribute. This is the case for the vector 
current, because the improvement term does not affect the conserved charge. 

Since the improvement of spectral quantities only depends on the improvement of the 
action, the effect of the vector current improvement is expected to be less signihcant for the 
effective mass at moderate distances, where the vector correlator is dominated by the rho 
meson. This has implications for e.g. the scale-setting procedure dehned through the effec¬ 
tive mass of the vector correlator ra. We explicitly observed the effect of the improvement 
in this quantity to be small. 
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